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ABSTRACT 

We explore the possibility of a local origin for ultra high energy cosmic rays (UHECRs). 
Using the catalogue of Karachentsev et al. including nearby galaxies with distances less 
than lOMpc (Local Volume), we search for a correlation with the sample of UHECR 
events released so far by the Pierre Auger collaboration. The counterpart sample se- 
lection is performed with variable distance and luminosity cuts which extract the most 
likely sources in the catalogue. The probability of chance correlation after penalizing 
for scans is 0.96%, which corresponds to a correlation signal of 2.6cr. We find that the 
parameters that maximize the signal are = 3.0°, I?max = 4Mpc and Mb = — 15 for 
the maximum angular separation between cosmic rays and galaxy sources, maximum 
distance to the source, and sources brighter than B-band absolute magnitude respec- 
tively. This implies a preference for the UHECRs arrival directions to be correlated 
with the nearest and most luminous galaxies in the Local Volume, while the angular 
distance between the cosmic ray events and their possible sources is similar to that 
found by The Pierre Auger Collaboration using active galactic nuclei (AGNs) within 
70-lOOMpc instead of local galaxies. We note that nearby galaxies with D < lOMpc 
show a similar correlation with UHECRs as compared to well-known particle accel- 
erators such as AGNs, although less than 20% of cosmic ray events are correlated to 
a source in our study. However, the observational evidence for mixed composition in 
the high-energy end of the cosmic ray spectrum supports the possibility of a local ori- 
gin for UHECRs, as CNO nuclei can travel only few Mpc without strong attenuation 
by the GZK effect, whereas the observed suppression in the energy spectrum would 
require more distant sources in the case of pure proton composition interacting with 
the CMB. 
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1 INTRODUCTION 

Ultra-high energy cosmic rays (UHECRs) may become a 
complementary probe of some astrophysical objects in ad- 
dition to observations in multiple wavelengths. In fact, the 
detection of these particles with enough statistics would rep- 
resent the awake of the development of multi-messenger as- 
trophysics. The Pierre Auger obse rvatory in Malargiie, Ar- 
gentina (see [Abraham et al.l |2004| ) has devoted a large ef- 
fort in the construction of an unprecedented array of water 
Cerenkov tanks covering an area of ~3,000 km^, together 
with four fluorescence telescopes which allow for increased 
accuracy in energy measurements. This effort has proven to 
be fruitful as the first relevant results were obtained even be- 
fore the completion of the entire experiment array. In partic- 



* E-mail: ajcv@iaa.es 



ular, the unprecedented statistics on the detection of ultra- 
high energy cosmic rays above 10EeV=10^^eU ~ IJ with 
arrival directions measured with an accuracy better than 1° 
has allowed to search for possible astrophysical sources of 
these parti cles with higher reliability than previous experi- 
ments fe.g. lFinlev fc Westerhogi2004l ). 

By the end of 2007, the Pierre Auger Collabora- 
tion concluded that the arrival directions of UHECRs 
above 55EeV correlate with the positions of nearby ac- 
tive galactic nuclei (AGN), or other objects which trace, 
in the same way as AGNs, the Large Scale Struc- 
ture of the Universe (The Pierre Auger Collabo rationll2007l . 
iThe Pierre Auger C ollaboration 2008). Statistics of these 
cosmic ray events has doubled since then (from 27 to 
58 events), and yet the correlation with these objects 
has not strengthened as compared to p revious estimations 
iThe Pierre Auger Collaboration! l2009bh . This supports the 



2 A. J. Cuesta & F. Prada 



fact that this scenario is not the only one which can re- 
produce the observed data, although increasing statistics of 
UHECR events may be able to finally disentangle their real 
origin. Expectations in the cosmic-ray com munity are opti- 
mistic, and as stated in lOlinto et alj (|2009l ) . this issue may 
find a solution in the next decade. 

Indeed, there has been a plethora of studies in the lit- 
erature trying to find out a correlation of the UHECRs with 
other candidate sources such as particular objects like galaxy 
clusters, GR Bs or other potential ho s ts of energetic phe- 
nomena (e.g. IPierpaoli fc Farrad |2005| . iMurase et al] |2008| . 
iGhisellini et al . 2008) and structures like the Supergalac- 
tic Plane (|Stanevl 1200^ . Some of these works did report 
a similar significance as compared to the Auger result on 
the correlation with the angular positions of nearby AGNs. 
These similar values of the probability of chance correlation 
failed to discriminate with enough robustness between dif- 
ferent cosmic-ray acceleration sites. In order to get stronger 
conclusions from this scarce set of UHECR events, there 
have been several efforts on the development of new meth- 
ods to estimate the chance probability of c orrelation in the 
case of isotropic flux with potential sou rces (|Ave et al.ll2009l . 
iThe Pierre Auger Collaborationll2009bl ). Despite all these ef- 
forts, there is no preferred candidate source yet. 

Therefore, any complementary information from other 
measurements, like the shape of the energy spectrum, 
or the composition of these cosmic rays, plays now a 
decisive role. For example, interactions wit h CMB pho- 
tons via the Greisen-Zatsepin-Kuzmin effect (|Greisenlll96^ . 
IZatsepin fc Kuz'minI Il966f ) are expected to decimate the 
statistics of the UHECRs arriving to the Earth as a function 
of their compositi on and travel distance (e.g. ICroninlboOSl . 
lOlinto et al.ll2009l). Recent claims of a detection o f this GZK 
cutoff d Abbasi et al]|2008l . [Abraham et"al]|2008l 'l raises the 
question of how far away these cosmic accelerators must be 
located in order to present the observed suppression. Be- 
sides, simple considerations regarding the coincidence of the 
GZK and UHECR clustering energies point to the fact that 
the composition of UHECRs must deviate from pure proton 
composition and tend to be CNO-like (jPermcr 2008), which 
is favoured by observations llBellido et al.ll2009i . lUnger et al.l 
I2OO9I . iThe Pierre Auger Collaboration! l2009al ). As CNO nu- 
clei can travel shorter distances as compared to protons, this 
sets a strong limit on the distance of UHECR sources, with- 
out suffering from the large deflection angl es by the galac- 
tic m agnetic held in the case of iron nuclei (|Takami fc Satol 
I2OO9I ). 

In this paper we address the possibility of a local origin 
for the Auger UHECRs. We use the on l y data published so 
far in iThe Pierre Auger Collaboration! l|2008l ) , which com- 
prises 27 UHECR events detected prior to 1 September 
2007, with energies above 57EeV (which has been revised 
to 55EeV due to new calibration procedures). This en- 
ergy scale turns out to be however ver y interesting in 
the se arch for these comic ray sources, as iMollerach et al.l 
l|2009l ) showed that the clustering above this energy is sta- 
tistically signiflcant. Most of the rece nt studies, including 
iThe Pierre Auger CoUaborationI (|2008l ). have focused so far 
on sources located eis far as 75 — lOOMpc, which corre- 
sponds to the GZK horizon for protons. Here we assess if 
there is any correlation signal with nearby objects with dis- 
tance D < lOMpc (Local Volume) which is roughly the GZK 



horizon for CNO nuclei. In order to explore this hypothesis, 
w e take advantage of t h e Loc al Volume catalogue of galaxies 
of lKarachentsev etal] (|2004| ). which is about 70-80% com- 
plete up to SMpc. This allows us to evaluate the likelihood 
of galaxies in the Local Volume as UHECR accelerator sites. 

This paper is organized as follows: in Section [5] we de- 
scribe the methodology used in our analysis and the cata- 
logue of nearby galaxies in the Local Volume. In Section [3] 
we present the results of the chance probability of the corre- 
lation of UHECRs with these objects. We discuss our results 
and conclude in Section |4] 



2 METHODOLOGY 

In this paper w e compare the distribution of UHECR arrival 
directions from lThe Pierre Auger Collaboration (2008) with 
the positi ons on the sky of nearby galaxies from the cat- 
alogue of lKarachentsev et al.l (|2004 ). This catalogue com- 
prises 451 neighbouring galaxies up to an estimated distance 
of D < lOMpc or a radial velocity VLg < 550km s~^ with 
respect to the Local Group centroid. In Figure [1] we rep- 
resent an Aitoff projection of the angular distribution of 
cosmic-ray arrival directions (circles) t ogether with the lo- 
cation s of all the galaxies in the whole iKarachentsev et al.l 
(j2004l ) Local Volume catalogue (stars), both in galactic coor- 
dinates. The Supergalactic plane is also shown with dashed 
line for reference. As expected, a large number of nearby 
galaxies in the catalogue are clustered near the Supergalac- 
tic plane. This makes these particular sources good can- 
didates for cosmic ray accelerators, as the arrival direc- 
tion s of UHECR s seem to be correlated with this plane 
(see IStanevll2008l ). In particular, there are a large number 
of these UHECR events which are clustered near the Su- 
pergalactic plane at the location of Centaurus A, as previ- 
ously pointed out bvlThe Pierre Auger CoUaborationI (|2008|) 
and conflrmed in iThe Pierre Auger Collaboration ( 2009lJ ). 
although some studies conclude that t his region is unlikel y 
to be a significant source of UHECRs l|Croston et al.ll2009l ). 

The informat i on a vailable in the catalogue of 
IKarachentsev et al] l|2004) includes the absolute magnitude 
in the _B-band and real distance information of the galaxies 
therein. This allows us to select the nearest and most lu- 
minous galaxies, which might be associated to the highest 
UHECR flux provided that there is a correlation between B- 
band luminosity and cosmic ray acceleration power. Hence, 
the relevant quantities we choose to conform our param- 
eter space are: the angular separation 1/) between cosmic- 
ray events and galaxy sources, the threshold in the distance 
to the candidate sources Dmax, and the maximum absolute 
magnitude in _B-band AIb- We did not study any depen- 
dence in energy of cosmic-ray events, as the data released 
by the Auger collaboration in 2007 only included 27 events, 
which is too sparse to make sub-samples with different ener- 
gies. We show in Table [T] the values in the parameter space 
explored in our study for these three parameters tp, Dmax 
and Mb- We note that the most finely scanned quantity is 
the angular separation tp, which is related to the deflection 
angle by the intervening mag netic fields. The sma l l num ber 
of galaxies in the catalogue of lKarachentsev et al. 1 1|2004 ) do 
not allow to perform a fine slicing in all the three parame- 
ters simultaneously, so we choose Aip = 0°.l for comparison 
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Figure 1. The spatial distribution in galactic coordinates of 
cosmic-ray arriv al directions (circles) and nearby galaxies from 
the catalogue of iKarachentsev et al.|[2004l (stars) shown in Aitoff 
projection. The Supergalactic plane is also shown with dashed 
line for reference. 



Table 1. Parameter space explored in this study. We show the 
minimum and maximum values, Jfmin and Xmax, as well as the 
scanning step AX, for the three parameters explored in our study, 
i.e., 1/1, Dmax and AfB(see text for details). 



Quantity X 1/) (°) D^ax (Mpc) Mg, 



-^min 


1.0 


3.0 


-18 




8.0 


10.0 


-6 


AX 


0.1 


0.5 


3 



with Fig. 3 in Ixhe Pierre Auger Collaboratioij l|2008l ). and 
took several slices in distance and luminosity. The sampled 
values according to Table [T] extends over a grid of only 5,325 
points in this space, which is computationally inexpensive. 

In order to have an estimation of the probability of re- 
jecting the null hypothesis (i.e. that the UHECRs in our 
sample show a distribution compatible with isotropy), we 
need to calculate the proba bility of correlation with a possi- 
ble iso tropic flux, following iThe Pierre Auger Collaboration! 
(|2008l ) . Provided that a particular combination of values for 
our three parameters can show a pronounced minimum in 
this probability, we have to explore the full parameter space 
so that we can study the case in which isotropy seems to be 
unlikely. Therefore, for each set of values of our parameters 
we evaluate the probability of isotropic flux associated to (at 
least) A^corr correlated source-cosmic ray event pairs. This 
is given as expected by the binomial distribution: 

Pdata= E (^.")p^(l-p)'"-- (1) 

where p is the probability that an event drawn from 
isotropy correlates with a galaxy within tp, using the sub- 
sample given by the other two parameters. To make this 
more clear, let us assume that the observatory has uniform 
exposure for the entire sky. In this case the circles in the 
sky of radius V' around galaxies in the catalogue which are 
brighter than Mb and nearer than Dmax deflne an area which 



is a fraction p of the full 4n sr sky. The value of p is esti- 
mated by MonteCarlo sampling of the sphere, calculating 
the fraction of points whose distance to a galaxy in the cat- 
alogue is smaller than a given angle ip. For very small values 
of ip, this estimation is straightforward as the areas covered 
by individual galaxies do not overlap, and hence the total 
area is equal to the number of galaxies A'^ times the area 
covered by the spherical cap subtended by the angle tp, i.e. 
0.5A''(1 — cos ■)/')■ On the other hand, for large values of 1/) or 
large number of galaxies, the area around individual galax- 
ies will overlap and this approximation is no longer valid, 
so a numerical estimation is needed. In the general case, we 
also need to take into account the declination dependence on 
the exposure due to partial sky coverage by the observatory. 
This is done by generating events according to the relative 
exposure distribution uj{5). The number of points we used 
in the MonteCarlo sampling of the sphere is 10*' so we have 
a 10"'^ uncertainty on the determination of p for each node 
in the grid of the parameter space. 



3 RESULTS 

We now apply the method described in the previous Section 
to calculate the probability in the case of isotropic flux of 
Acorr correlated source-cosmic ray event pairs, as a func- 
tion of the parameters tp, Dmax and Mb- A full exploratory 
scan is performed over the entire parameter space, and as 
a result we find that this probability reaches its absolute 
minimum (hereafter Pdata) at cosmic ray event deflections 
below tp = 3.0° and galaxies accelerating UHECRs with 
distances below Dmax = 4Mpc and brighter than Mb=- 
15. This implies a preference for the UHECRs arrival di- 
rections to be correlated with the nearest and most lumi- 
nous galaxies in the catalogue, while the angular distance 
between the cosmic ray e vents and their possible sources 
is similar to that found bv lThe Pierre Auger CollaboratiorJ 
(2008) using AGNs instead of local galaxies. In order to 
visualize the dependence of this probability with individ- 
ual parameters, we refer to Figure (2] Here we display the 
chance probability of correlation with the Auger UHECRs 
above E > 55EeV as a function of each one of the scan 
parameters, keeping the other two parameters fixed at the 
absolute minimum of the probability of isotropic fiux. When 
the angular distance ^ varies, we find several local min- 
ima, although the region in which the chance correlation 
probability is below 0.1% is in the range of 2-4°. On the 
other hand, the behaviour with the other two scan param- 
eters is quite robust with a single minimum, although this 
is also due to the coarseness in the scanning values in this 
case. Nevertheless, it seems clear that including distant and 
intrinsically faint sources decreases the correlation signal. 
Moreover, given the value of tp the deflection of cosmic ray 
trajectories due to intervening magnetic fields is not large, 
suggesting that these events are hgh t nuclei, as suggested by 
iThe Pierre Auger CoUaborationI (12008). 

As shown in Figure [21 the probability at the absolute 
minimum is Pdata ~ 1-1 x 10'*^. This signal is simila r to th e 
previous result of iThe Pierre Auger Collaboration! (|2008l ) , 
who found Pdata = 2 X 10""^, and even stronger t han th e 
revised result by iThe Pierre Auger Collaboration! (|2009bl ) , 
who found that Pdata = 6 x 10~^ for the correlation with 
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Figure 2. The chance correlation probabiUty as a function of each one of the scan parameters, holding the other two fixed at their 
minimum of Pdata- ^^fl panel: Angular separation ■(/>. Middle panel: Maximum galaxy distance Dmax- Right panel: Maximum absolute 
magnitude in the _B-band. 



AGNs up to 75Mpc. Nevertheless, our result is not free from 
the negative impact of trial factors in a posteriori anisotropy 
searches, as it was not tested on an independent data set 
with the scan parameters specified a priori. For reference, the 
Auger result including all the 27 events is Pdata ~ 4.6 x 10~^. 
Also, our value is very sensitive to the low number of cor- 
relating events at this minimum (five cosmic rays associ- 
ated to four sources), so removing one UHECR can boost 
^tiata to 1.10 X 10~^. However, it is important to note that 
Eq. m is not the true probability o f chance correlation 
due to isotropic fiux, as pointed out bv lFinlev fc Westerhofj 
(|2004[) . Although this exploratory scan is useful in order to 
detect the location of the strongest potential correlation sig- 
nal, the fact that we have chosen our parameters in such a 
way to get the minimum value of this probability has to 
be taken into account. In order to correct the value Pdata 
that we have just obtained, we need to include a proba- 
bility penalty due to parameter selection. This is fixed by 
performing umc MonteCarlo realizations of the cosmic-ray 
data set, which are drawn from a isotropic distribution but 
they are also selected in such a way that they mimic the 
events which the Pierre Auger Observatory might detect, 
due to declination dependence of the exposure. Therefore, 
the mock cosmic ray data set assume isotropic flux but the 
probability of a given arrival direction (a, 5) is proportional 
to the exposure function of the Au ger observatory on that 
position (see e.g. ICuoco et al.|[2006l ). 

The isotropic samples are obtained as usual by generat- 
ing random points on the sphere, i.e. picking uniformly dis- 
tributed random numbers in sin 5 and a, so that the number 
density of points is completely homogeneous over the solid 
angle unit dfJ = cos 5A5Aa. To account for the Auger ob- 
served sky, we accept a random point on the sphere if an 
uniformly distributed random number between and 1 is 
smaller than the relative exposure u normalized by its max- 
i mum value c Jmax, o therwise this point is rejected. According 
to lSommersI (|200ll ). the relative exposure is given by: 

uj{5) oc cos(A) sin(Qm(C)) cos((5) + Q.m{C,) sin(A) sin((5) (2) 

where A ~ —35° is the latitude corresponding to the location 



of the Pierre Auger Observatory, C, is given by: 
cos(0max) — sin(A) sin((5) 



cos(A) cos((5) 



(3) 



(^max = 7r/3 is the cut in zenith angle 6 of the arrival direc- 
tions applied by Auger), and a™ is defined as follows: 



Q-m = { arccosi^ 




if C < ~i 

if - 1< C < +1 
if C > +1 



(4) 



Therefore, if we perform the same exploratory scan in 
the same way as in the beginning of this Section, for each 
different set i of Monte Carlo generated arrival directions 
{i = 1, . . . , nMc), we can just count the number of trials njjc 
for which the absolute minimum (over the entire parameter 
space) of the probability of isotropic fiux is lower than the 
actually observed minimum, - P^in ^ fdata- The chance p rob- 



! p r 

abfiity of observing Pdata is (|Finlev fc WesterhofjlioM ): 

P=^ (5) 

nyic 

We can consider this value as an estimate of the 
chance probability of correlation due to isotropic fiux, as in 
iThe Pierre Auger CoUaborationI |200i). We performed 10 
Monte Carlo realizations of the Auger cosmic ray data set, 
in order to keep our uncertainties below 0.1%. The value 
we obtain for this probability after correcting for parameter 
selection is the following one: 



Pchancc = (9.6 ± 1.0) X 10" 



(6) 



We find that the probability of chance correlation of 
the highest energy (E > 55EeV) Auge r events with nearby 
galaxies from the iKarachentsev et al.l l|2004h catalogue is 
around 0.96%. This small probability corresponds to a 2. Go- 
correlation between both samples. We point out that more 
than the half of the 27 Auger events above E > 55EeV are 
within 4° of a nearby {D < lOMpc) galaxy, which made it 
worth of exploring. This result goes along the same direction 
as the previous studies, finding a hint of violation of isotropy 
in the Auger sample. 

There are 5 out of 27 correlating events, while only 0.5 
should be expected on average in the case of isotropic fiux. 
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as the effective fraction of the sky for the values of the pa- 
rameters at the minimum is p = 0.018. Of course the number 
of correlating events is much low er than the value found by 
iThe Pierre Auger Collaboration! (2008) (20 out of 27, with 
5.6 expected in the case of isotropic flux), but it should be 
considered that the effective fraction of the sky covered (at 
the minimum Pdata) by 442 objects in the Veron-Cetty AGN 
catalogue with z ^ 0.017 is much larger (p ~ 0.21) than the 
31 objects with our parameters (p = 0.018). 

The galaxies which correlate with the arrival directions 
of the cosmic rays, are NGC300 (spiral galaxy type SA), 
NGC4945 (Sy2), NGC5102 (HII) and NGC5128 (Sy2). With 
exception of NGC300, they mostly show nuclear activity and 
both S eyfert galaxies are included in t h e AG N catalogue 
used bv lThe Pierre Auger Collaboration! (l2008t l. In particu- 
lar NGC5128, which is associated with 2 cosmic-ray events, 
is inside the well known Cen A region. Thus, it might hap- 
pen that the correlation we find could just be inherited by 
the correlation found by the Auger Collaboration if the lat- 
ter is finally confirmed. This is reasonable, as only a few 
tens of AGNs from the Veron-Cetty catalogue are located 
at a distance of D < lOMpc, but it is important to keep 
in mind that CNO UHECRs from larger distances can be 
heavily suppressed. 



4 DISCUSSION AND CONCLUSION 

We have explored the possibility of the correlation of nearby 
galaxies from the catalogue of JCarachontsev ct al. (2004) 
with ultra high energy cosmic ray events from the Pierre 
Auger Collaboration. This analysis is highly motivated in 
the light that a mixed CNO (as opposed to pure proton) 
composition of UHECRs, which is the preferred scenario in- 
dicated by current data, could not survive the GZK cutoff 
unless the sources are located in the Local Volume. There- 
fore, any correlation with these nearby sources would ex- 
plain both composition and correlation at the same time. 
The analysis of the data concludes a 2.6(t correlation be- 
tween both samples. There are 5 out of 27 correlating events, 
while 0.5 should be expected on average in the case of 
isotropic flux. Due to the small number of correlating events, 
we cannot specify the scan parameters a priori, thus mak- 
ing our results sensitive to the negative impact of a pos- 
teriori anisotropy searches. The probability of chance cor- 
relation due to isotropic flux after penalizing for scans is 
^chance = 0.96 X lO"'^, rejecting the hypothesis of isotropy 
of arrival directions of UHECRs at 99% confidence level. 
For comparison, the same result before penalizing for scans 
is ~ 1 X 10~* which is similar to the ~ 2 x 10~* result 
(with sca n parameters set a priori) from the c orrelation 
found by iThe Pierre Auger Collaboration! (120081 ) with the 
AGNs in the Veron-Cetty catalogue, and stronger than the 
revised value of 6 x 10~^ when they include the new 
UHECR events from 1 September 2007 to 31 March 2009 
ijThe Pierre Auger CoUaboration! 12009^ 1. This is an indi- 
cation that there is no clear source of the UHECRs, and 
finding a lower value of the probability of chance corre- 
lation does not guarantee that this is actually the pre- 
ferred source as compared to other sources which present 
similar values of Pchance . This was already noticed even 
in early papers such as iTinvakov fc Tkachev! |20o3), and 



in any case it seems un likely to get a > 3.3(t detection 
(jKashti fc Waxman 20081). Interestingly , most of the objects 
in the catalogue of iKarachentsev et al. I (Hooi) which corre- 
late with Auger events, are already included in the AGN 
catalogue, which has around 30 objects with D < lOMpc. 
This may indicate that the relation between nearby galaxies 
and UHECRs found in our work c ould be inherited from the 
iThe Pierre Auger C ollaboration' (2008') result. In any case, 
the so-called AGN hypothesis has proven to be c ontrover- 
sial ijCorbunov et al ] |2008l . lMoskalenko et al ] |2009l ). stressed 
by the deficit of events from the direction of the Virgo su- 
percluster. On the other hand, the correlation with the dis- 
tribution of luminous ma tter seems to be well established 
(jKashti fc Waxman!l200d ) , without making any reference to 
any particular source which can accelerate particles up to 
these energies. 

The small number of correlating events with nearby 
galaxies might also point to the possibility of the absence 
of UHECR accelerators in our local environment, although 
there is no de finitive answer yet. Surpris ingly, most of the 
objects in the IKarachentsev et al.l (|2004h catalogue which 
correlate with the UHECR sample (with the exception 
of NGC300) are located in a small region in the sky, 
with Supergalactic coordinates 150° < SGL < 165° and 
— 10° < SGB < —5°. This region is near the Supergalactic 
plane, pointing to preferential arrival directions from nearby 
galaxies located in a high density region. Future measure- 
ments and hopefully the next data release with a total of 
58 events as announced in lThe Pierre Auger Collaboration! 
(2009b) will help to clarify the origin of nearby and distant 
UHECRs, which will allow to study more in detail the 
acceleration physics behind them. 
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